Atomic parity nonconservation in Ra+ 
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We report on a theoretical analysis of the suitability of the 7s^Si/2 ^^ &d^D^/2 transition in 
singly ionized radium to measure parity nonconservation, in the light of an experiment planned 
at the KVI of the University of Groningen. Relativistic coupled-cluster theory has been employed 
to perform an ab initio calculation of the parity nonconserving electric dipole amplitude of this 
transition, including single, double, and leading triple excitations. We discuss the prospects for a 
sub-1% precision test of the electroweak theory of particle physics. 
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In atomic systems, parity is broken due to the exchange 
of the neutral vector boson Z'^ , which mediates the weak 
interaction between the atomic electrons and the quarks 
in the nucleus. This atomic parity nonconservation 
(APNC) gives rise to small parity nonconserving electric 
dipole transition amplitudes (-BlpNc)- The APNC effect 
gets strongly enhanced in heavy atoms and can be mea- 
sured by the interference of i?lpNC with a suppressed 
electromagnetic transition amplitude {Ml, E2) [H [2]. 
The accurate measurement of the 63^5*1/2 <-^ 7s ^5*1/2 
transition in atomic ^'^'^Cs by the Boulder group [31 H] 
constitutes a precision test of the electroweak sector of 
the Standard Model (SM) of particle physics [5]. By 
combining the measurement with a many-body atomic 
structure calculation, the weak nuclear charge could be 
determined [6\. 

The importance of APNC to particle physics is a 
strong incentive to further pursue these challenging ex- 
periments. With the experimental and theoretical accu- 
racies at an impressive 0.35% |3l H] and 0.5% ^ level, 
respectively, the ^'^^Cs result agrees with the SM predic- 
tion within one standard deviation. Nevertheless, it is 
desirable to consider other candidates for APNC stud- 
ies, see e.g. Ref. [7]. New experiments have been pro- 
posed for Cs J8| and Fr [9j atoms. Of special interest is 
the proposal by Fortson to measure APNC in one single 
laser-cooled and trapped ion [5j. Such single-ion exper- 
iments offer important benefits, such as long coherence 
times and precise control of various systematic effects. 
Promising ions from the experimental and atomic-theory 
point of view are heavy alkali- like ions, in particular Ba'^ 
and Ra"*" |10j . Proof-of-principle experiments have been 
carried out with ^^^Ba+ by the Fortson group [TUIITTIIT^ . 
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FIG. 1: Relevant energy levels in Ra^ 



At the TRI/iP facility ^^ [TJ at the accelerator in- 
stitute KVI in Groningen an APNC experiment on Ra"*" 
is in progress [15]. An important advantage of Ra"*" is 
that all relevant transitions are in the optical regime, c/. 
Fig. 1, and thus are accessible by commercially available 
solid-state laser technology. The goal is to measure the 
-ElpNC amplitude of the ls'^Si/2 <-> &d^D^/2 transition. 
We address here the question what the prospects are to 
push the corresponding atomic theory below 1%, such 
that the experiment can serve as a high-precision test of 
the SM. We analyze various relevant properties of Ra"*" 
and assess the remaining uncertainties. 

The parity-nonconserving nuclear-spin independent 
(NSI) interaction is due to the electron-quark neutral 
weak interaction, the Hamiltonian of which is given by 
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where Gp is the Fermi constant, £»nuc the nuclear density. 



and 75 is the standard Dirac matrix; Qw is the weak nu- 
clear charge, which is equal to {2Z + N)ciu + {"^N + Z)cid 
in terms of the coupling constants of the electron to the 
up and down quarks; Z and N are the number of protons 
and neutrons. The Hamiltonian in Eq. ([l]) mixes atomic 
states of opposite parity but with the same angular mo- 
mentum. Its strength is weak enough to consider it as a 
first-order perturbation. We therefore write the valence 
state {v) atomic wave function as 



|vI/.) = |vI/(o)) + G^|vI/W) 
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amphtude between the ls^Si/2 (= i) and 6(i D3/2 {— /) 
states is given by 
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which, after expansion, takes the form 
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where D is the electric dipole {El) operator, I and J 
represent the allowed intermediate states, and E is the 
energy of the state. An accurate determination of i^lpNc 
depends on the precision of the matrix elements of D 
and of Hp^\., and of the energy differences between the 
different states. At the same time, it is also important 
to take all intermediate states into account, something 
which is not possible in the often-used sum-over-states 
approach. We therefore employ the relativistic coupled- 
cluster (RCC) theory, which allows us to evaluate the 
properties to all orders in perturbation theory. The 
RCC method was previously used to calculate APNC in 
"7Ba+ with sub-1% accuracy [16] . 

We obtain the first-order wave functions of Eq. (Is]) in 
the RCC framework as the solution of 

GpiHuc " E,)\^i'^) = -i/^^^cl^f > , (5) 

where v stands for valence electron, which is either i or 
/. The unperturbed and perturbed wave functions are 
expressed as 



|vI/W)=exp(T(")){l-H5(0)}|a>,), 



(6) 



and 



|vI/W) = eMT'){l + S',}\^,) 

= exp(T(o))(T(i){l-HS(°)} + {5W})|a>.) ,(7) 



respectively, where |$„) is the mean-field wave function 
determined with the Dirac-Fock (DF) method. T and 
Sv are the core and valence-core RCC correlation oper- 
ators, respectively, where the superscript indicates in 
the presence of the Coulomb interaction, the prime (') in- 
dicates in the presence of both the Coulomb and APNC 
interaction, and 1 indicates their linear approximations. 
Substituting the above expressions in Eq. pi, we obtain 
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where jvpi, ) is the atomic wave function of the Dirac- 
Coulomb (DC) Hamiltonian (i^Dc) and \'9v ) is the first- where 
order correction due to the PNC NSI interaction. 

To a first-order approximation, the _ElpNc transition 

, 20 , /_ „\ „„J R^2r 
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and Di°) = exp{T'-°'^^)Dexp{T'-°'^). The matrix element 
is evaluated using the generalized Wick's theorem [T5] . 

Our RCC work has two salient features. We evaluate 
-El PNC by using the direct solution of the first-order per- 
turbed equation as given in Eq. ([5]) rather than summing 
over a finite number of intermediate states [T7] . The core 
correlation effects modified by the parity nonconserving 
weak interaction are evaluated to all orders through T^^' 
in the framework of the relativistic CCSD(T) method. 
These effects cancel strongly in Cs and Fr, where both 
the initial and final states are 5-states. However, it is 
essential to consider them accurately in the S-D tran- 
sitions in Ba+ and Ra+, where these contributions are 
significant. For our calculation, we have used numerical 
DF/V^~^ orbitals to describe the occupied and bound 
virtual orbitals. The continuum states were represented 
by V^^^ Gaussian-type orbitals (GTOs) [18 using the 
parameters a = 5.25 x 10~^ and (3 = 2.73. The finite 
size of the nucleus is accounted for by assuming a Fermi 
charge distribution [T5]. 

In Table |l] we present our RCC results for the £^lpNC 
amplitude of the 7s '^Si^ <-> 6d ^0^/2 transition in the iso- 
tope ^^^Ra+. Shown are the results of the DF method, 
of the RCC method with single and double excita- 
tions (CCSD), and with the leading triple excitations 
(CCSD(T)). The difference between the CCSD(T) and 
CCSD results is small. Our best value is the CCSD(T) 
result -ElpNC = 46.4 x lO'^^ieao{-Qw/N). Also shown 
are two results of Dzuba et al. [12] : in a sum-over-states 



TABLE I: BlpNC for the 7s ^Si/2 <-^ 6d ^L>3/2 transition in the 
isotope ■^■^^Ra^, in units of W^^^ieao(—Qw/N). 



This work 



Ref. nn 



DF 40.4 

CCSD 46.1 Mixed-states 42.9 

CCSD(T) 46.4 Sum-over-states 45.9 



approach they found 45.9 x 10~^^ieao{—Qw /N)', in a 
mixed-states approach, wherein the APNC interaction 
explicitly mixes states of opposite parity, they obtained 
42.9 X 10~^^ieao{—Q\Y/N)- Neither calculation includes 
structural radiation, the weak correlation potential, and 
normalization of states, effects which are included by us. 

The E'lpNC amplitude for the 6s ^<S'i/2 <-> 7s ^51/2 tran- 
sition in Cs is about 0.9xl0~^^?eao(— Qiy/iV) Hj. Thus, 
the APNC effect in Ra+ is larger by a factor close to 
50. In heavy atoms, ANPC gets enhanced by the overall 
factor KrZ'^QwiZ,N), where Qw ^ N ^ Z, and Kr is 
a relativistic factor that depends on the nuclear charge 
and radius. This is the "faster-than-Z'^ law" [T], which 
implies that Ra+ is favored over Cs by a factor of about 
20. An additional factor of around 2 can be understood 
as follows. For Cs (and Fr) the S-S transition is used, for 
Ra+ (and Ba"'") the S-D transition. Since the Z^-hoson 
is very heavy, the weak interaction between electrons and 
the quarks in the nucleus has (almost) zero range. The 
overlap of the electrons with the nucleus is largest for the 
5' states, and thus the mixing of the P-states into the S 
states gives the major contribution to iJlpNC- However, 
in Cs and Fr the initial and final S states contribute with 
opposite signs, which leads to a significant cancellation 
in SlpNC- In fact, for Cs there are three dominant terms 
in the sum over the states, which add up to a total value 
that is half the size of the largest individual term [^ . 

The S-D transitions in Ba"*" and Ra+ do not suffer 
from such a cancellation, since the contribution from 
the D-state to APNC is small. In Tabic |TT] wc ana- 
lyze which intermediate states contribute most to the 
total sum. Clearly, in contrast to the Cs S-S case, 
the sum is strongly dominated by one term: the con- 
tribution from the 7p^Pif2 state. These qualitative re- 
sults are robust, they are consistent with the findings 
of Ref. [in], and they are, in fact, already borne out by 
a simple calculation with quantum-defect theory, analo- 
gous to Ref. HO]. This simple estimate gives for Cs, Fr, 
Ba"*", and Ra"*" results accurate to some 10%; for Ra+ we 
find £^lpNC = 45(4) x lO~^Heaoi-Qw/N). 

In Table |III| we present our results for the excitation 
energies. El transition amplitudes, and hyperfine con- 
stants for the relevant transitions and states in Ra"*". We 
also list experimental values where available. For the ex- 
citation energies, we compare to the only available spec- 



TABLE IT: The contributions to the _BlpNC from the different 
P-states (%). 
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up J^3/2 

7p^P\/2 bound 
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8.7 
-15 

111 
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bound 

bound 

continuum 

continuum 



-3.3 

-0.7 

-0.1 

1.1 



TABLE III: Excitation energies, El transition amplitudes, 
and Ai /gi for different low- lying states of Ra+. 



Transition 


7s^5i/2 

7p'Pl/2 


73^51/2 
7p'P3/2 


6d^D3/2 

7p'Pl/2 


6d^D3/2 
7p'P3/2 


Excitation energy [cm"""^] 
This work 21509 
Experiment [2lj 21351 


26440 
26209 


9734 
9267 


14665 
14125 


El transition amplitude [a.u.] 
This work 3.31 
Ref. |19| 3.223 
GTOs [22] 3.28 


4.58 

4.477 

4.54 


3.68 

3.363 

3.64 


1.56 

1.504 

1.54 


State 


7s ^5'i/2 


7p'Pl/2 


7p'P3/2 


6d2p)3/2 



Hyperfine interaction constant (Ai/gi) [MHz] 

This work 19689.37 3713.75 312.91 441.67 

Experiment [23] 18772 3691 314.12 



troscopy measurement [21], which dates back to 1933. 
For the El transition amplitudes, for which there are no 
experimental data, we list the results of Ref. [19] and of 
our previous work [23 using GTOs. Therein, the life- 
times of the metastable Z?-states were calculated to be 
0.627(4) s for Gd^D^/^ and 0.297(4) s for M^D^/2. 

Since the hyperfine structure is a good probe of the 
wave functions at the nucleus, we have, in order to esti- 
mate the accuracy of the H^^q matrix elements, cal- 
culated the ratio between the magnetic dipole hyper- 
fine structure {A) and the nuclear gyromagnetic {g) con- 
stants, neglecting isotope effects, and compared these 
with experimental results for Ra+ from ISOLDE [33] JH] . 
Our calculated value for [Ai/gi{7Si/2)Ai/ gi{7Pi/2)Y/'^ 
differs by 3% from the experimental value, which is 
a reasonable estimate for the dominant uncertainty in 
the atomic theory. Thus, our best value for the par- 
ity nonconserving El amplitude in Ra+ is _BlpNc ~ 
46.4(1.4) X 10-"ieao(-Qw/iV). 

It appears feasible to push the accuracy of the atomic 
theory for Ra"*" to the sub-1% level. Improvements along 
several lines are in progress. The Breit interaction [^[^ 
and QED corrections, which contribute around 1%, need 
to be included. The neutron-skin effect [27], which also 
contributes at the sub-1% level, has to be investigated. 
However, at the same time it is clear that experimental 
information to test the atomic theory is severely lacking. 
Not all relevant energy levels are known [5T], there is no 
experimental information on the El transition strengths, 
nor on the lifetimes of the D-states. It is highly desirable 
to have more experimental data on the magnetic dipole 
and electric quadrupole hyperfine constants A and B for 
the various Ra"*" isotopes. The extraction of these con- 
stants [231 123] is model dependent, and ideally one would 
like to use a single consistent ah initio framework for this. 



TABLE IV: The properties of the isotopes of Ra^ suitable for 
a single-ion APNC experiment. A is the mass number, I is 
the nuclear spin, and T1/2 the half-life time. 



computer in C-DAC, Bangalore. 



A 



Tl/2 



Possible production reaction 



-^ ^"Ra + 7n 

22^Ra + -^X + an + 6p 

^^''Ra +'^X+an+bp 



213 1/2" 2.74(6) min ^osp^ _^_ i2q 

223 3/2+ 11.43(5) d p + ^^^Th - 

224 0+ 3.6319(23) d p + ^^^Th - 

225 1/2+ 14.9(2) d ^^^Th -^ ^^^Ra + a 

226 0+ 1600(7) y Commercially available 

227 3/2+ 42.2(5) min p + ^^^Th -^ ^^^Ra +'^X + an + 6p 



At the TRI/iP facility, radium isotopes can be pro- 
duced in fusion and evaporation or spallation reactions. 
The ions can be collected in a radio- frequency trap where 
they can be laser-cooled on the 7s ^51/2 *-^ 7p ^Pi/2 reso- 
nance line at 468 nm, with repumping via the 6d ^1)3/2 ^^ 
7p^Pi/2 transition at 1.08 /xm, for which strong lasers are 
available. They will then be transferred to a miniature 
trap for the single-ion experiment, where techniques sim- 
ilar to Refs. [ini HH [H] will be applied to perform the 
measurements. In particular, the 7s'^Si/2 ^^ '!p^P^/2 and 
6^^1)5/2 ^^ Tp'^P^/2 transitions at 381 nm and 801 nm, 
respectively, can be used for "shelving" [TH] . 

A list of the Ra+ isotopes suitable for a single-ion ex- 
periment is shown in Table |IV[ A half-life of the order of 
seconds is required for a high-precision single-ion experi- 
ment, but, on the other hand, it should not exceed a few 
days, so as to avoid long-lived radioactive contamination 
of the core equipment. Good candidates therefore are the 
odd isotopes ^^^Ra+ and ^^^Ra+, and the even isotope 
^^''Ra+. (Since the odd isotopes have a nonzero nuclear 
spin, the nuclear spin-dependent weak interaction will 
contribute to ii^lpNC [28]). The isotopes listed can all be 
produced at TRI/iP, and the experimental data required 
to constrain the atomic theory can be measured there 
with laser spectroscopy. Since multiple Ra+ isotopes will 
be available, the possibility exists to measure APNC in a 
chain of isotopes, which can help to eliminate remaining 
uncertainties in the atomic theory |29j . 

In conclusion, Ra+ appears to be an excellent candi- 
date for an APNC experiment, since SlpNC is large, the 
required lasers are all at convenient wavelengths, and 
one can exploit the high-precision techniques of single- 
ion trapping. The atomic theory needed for the interpre- 
tation of the experiment could reach an accuracy better 
than 1%, but precise experimental data for the relevant 
atomic properties are mandatory to achieve such a bench- 
mark. The prospects for APNC in Ra+ as a precise test 
of the SM look promising. 
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